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ABSTRACT 

We directly construct model-independent mass profiles of galaxy clusters from combined weak-lensing 
distortion and magnification measurements within a Bayesian statistical framework, which allows for a 
full parameter-space extraction of the underlying signal. This method applies to the full range of radius 
outside the Einstein radius, and recovers the absolute mass normalization. We apply our method to 
deep Subaru imaging of five high-mass (> 10 15 Af Q ) clusters, A1689, A1703, A370, C10024+17, and 
RXJ1347-11, to obtain accurate profiles to beyond the virial radius (r v ir)- For each cluster the lens 
distortion and magnification data are shown to be consistent with each other, and the total signal- 
to-noise ratio of the combined measurements ranges from 13 to 24 per cluster. We form a model- 
independent mass profile from stacking the clusters, which is detected at 37<r out to R rss 1.7r V j r . The 
projected logarithmic slope 72d(-R) = dlnE/dlni? steepens from 72D = —1-01 ± 0.09 at R ss 0.1r v i r 
to 72D = —1-92 ± 0.51 at R « 0.9r V i r . We also derive for each cluster inner strong-lensing based 
mass profiles from deep Advanced Camera for Surveys observations with the Hubble Space Telescope, 
which we show overlap well with the outer Subaru-based profiles and together are well described by a 
generalized form of the Navarro-Frcnk- White profile, except for the ongoing merger RXJ1347-11, with 
modest variations in the central cusp slope (— din p/dhir ;$ 0.9). The improvement here from adding 
the magnification measurements is significant, ~ 30% in terms of cluster mass profile measurements, 
compared with the lensing distortion signal. 

Keywords: cosmology: observations — dark matter — galaxies: clusters: individual (A1689, A1703, 
A370, C10024+1654, RXJ1347-1145) — gravitational lensing: weak — gravitational lens- 
ing: strong 



1. INTRODUCTION 

Galaxy clusters provide an independent means of ex- 
amining any viable model of cosmic structure forma- 
tion through the growth of structure and by the form 
of their equilibrium mass profiles, complementing cos- 
mic microwave background and galaxy clustering obser- 
vations. A consistent framework of structure formation 
requires that most of the matter in the Universe is in the 
hitherto unknown form of dark matter, of an unknown 
nature, and that most of the energy filling the Universe 
today is in the form of a mysterious "dark energy" , char- 
acterized by a negative pressure. This model actually 
requires that the expansion rate of the Universe has re- 
cently changed sign and is currently accelerating. 

Clusters play a direct role in testing cosmological mod- 
els, providing several independent checks of any viable 
cosmology, including the current consensus A cold dark 
matter (ACDM) model. A spectacular example has 
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been recently provided from detailed lensing and X-ray 
observations of the "Bullet Cluster" (aka, IE0657-56; 
Markevitch ct al. 2004; Clowe ct al. 2006), which is a 
consequence of a high-speed collision between two clus- 
ter components with a mass ratio of the order of 6 : 1 
(Mastropictro & Burkcrt 2008), displaying a prominent 
bow shock preceding a cool bullet lying between the 
two clusters, implying these clusters passed through each 
other recently (Markevitch et al. 2002, 2004). Here the 
Bullet system reveals lensing mass contours that follow 
the bimodal distribution of cluster members, demon- 
strating that the bulk of the dark matter is relatively 
collisionless as galaxies (Clowe et al. 2006), as also shown 
by a comprehensive analysis of galaxy and dark-matter 
dynamics for A1689 (Lemze et al. 2011). Other cases of 
merging systems show that in general displacement of 
the hot gas relative to the dark matter is related to in- 
teraction (Jee et al. 2005; Okabe & Umetsu 2008). For 
dynamically-relaxed clusters, the form of the equilibrium 
mass profile reflects closely the distribution of dark mat- 
ter (see Mead et al. 2010) which, unlike galaxies, does 
not suffer from halo compression by adiabatic contraction 
of cooled gas. The majority of baryons in clusters are in 
the form of hot, diffuse X-ray emitting gas, and repre- 
sents only a minor fraction of the total lensing mass near 
the centers of clusters (Lemze et al. 2008; Umetsu et al. 
2009). 

The predicted Navarro-Frcnk- White profile (hereafter, 
NFW; Navarro ct al. 1996, 1997) derived from simula- 
tions based on collisionless, cold (non-relativistic) dark 
matter has a continuously-declining logarithmic gradi- 
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ent 73dM = dlnp/dhir towards the center of mass, 
much shallower than the isothermal case (730 = —2) 
within the characteristic scale radius, r s ( ;$ 300kpc/i _1 
for cluster-sized halos). A useful index of the degree 
of concentration, c v i r , compares the virial radius, r vu ., 
to r s of the NFW profile, c v j r = r y - 1T /r 8 . This has 
been confirmed thoroughly with higher resolution sim- 
ulations (Ghigna et al. 1998; Fukushige & Makino 1997; 
Okamoto & Habe 1999; Power et al. 2003; Navarro et al. 
2004; Neto et al. 2007), with some intrinsic variation 
related to the individual assembly history of a clus- 
ter (Jing & Suto 2000; Tasitsiomi et al. 2004). Grav- 
itational lensing observations are underway to pro- 
vide reliable and representative cluster mass profiles to 
test this since the first careful measurements showed 
that the NFW profile provides a good fit to the en- 
tire mass profile when weak and strong lensing are 
combined (Gavazzi et al. 2003; Broadhurst et al. 2005a, 
2008; Umetsu et al. 2010). Other well studied clusters 
with similarly high quality data are also in good agree- 
ment providing strong support for the CDM scenario 
(e.g., Okabe et al. 2010). 

Interestingly these studies reveal that although the 
dark matter is consistent with being cold, the predicted 
profile concentration of the standard ACDM model falls 
short of some lensing results (e.g., Broadhurst et al. 
2005a, 2008; Oguri et al. 2009). This observed tendency 
for higher proportion of mass to lie at smaller radius 
in projection is also indicated by the generally large Ein- 
stein radii determined from strong lensing of well studied 
clusters (Broadhurst & Barkana 2008; Zitrin et al. 2011) 
finding a substantial discrepancy with the predictions de- 
spite careful accounting for potential selection biases in- 
herent to lensing (Hennawi et al. 2007; Meneghetti et al. 
2010a). These observations could suggest either substan- 
tial mass projected along the line of sight, perhaps in 
part due to halo triaxiality (Oguri et al. 2005), or a large 
overconccntration of mass; the latter could imply modi- 
fication within the context of the CDM family of models. 

The abundance of massive clusters is very sensitive 
to the amplitude of the initial mass power spectrum 
(Shcth et al. 2001) representing the most massive objects 
to have collapsed under their own gravity, and confirmed 
by TV-body simulations of Hubble volumes (Evrard et al. 
2002). Such calculations predict for example that the 
single most massive cluster to be found in the universe 
is expected to be with Af v i r = 4 x 10 15 M Q out to 
z = 0.4 (see Figure 5 of Broadhurst & Barkana 2008), 
similar to the most massive known clusters detected lo- 
cally (Broadhurst et al. 2008). 8 At higher redshifts this 
comparison becomes more sensitive to the cosmological 
model, with an order of magnitude decline in the abun- 
dance of 10 15 M Q clusters at z > 0.8 compared to the 
present (Evrard et al. 2002). Hence, the existence of 
such massive clusters like XMMUJ2235-25 at z = 1.45 
(Jee et al. 2009), from lensing work, begins to motivate 
alternative ideas such as departures from Gaussian ini- 
tial density fluctuation spectrum, or higher levels of dark 
energy in the past (Sadeh & Rephaeli 2008), although 
some non-Gaussian models can be ruled out by using the 
cosmic X-ray background measurements (Lemze et al. 

8 A370 at z = 0.375 is currently the most massive known cluster 
measured reliably by lensing, M vir = (2.9 ± 0.3) X 10 15 M Q . 



2009b). 

The main attraction of gravitational lensing in the 
cluster regime (e.g., Hattori et al. 1999; Umetsu et al. 
1999; Bartelmann & Schneider 2001; Umetsu 2010) is the 
model-free determination of mass profiles allowed over a 
wide range of radius when the complementary effects of 
strong and weak lensing are combined (Broadhurst et al. 
2005a; Umetsu & Broadhurst 2008; Merten et al. 2009; 
Umetsu et al. 2010; Meneghetti et al. 2010b; Zitrin et al. 
2010). In practice, the quality of data required challenges 
with few facilities, which are able to generate data of suf- 
ficient precision to provide a significant detection of the 
weak lensing signal on an individual cluster basis. 

In this paper we aim to pursue in greater depth 
the utility of massive clusters for defining highest- 
precision mass profiles by combining all lensing in- 
formation available in the cluster regime. In par- 
ticular, we shall make a full use of the magnifica- 
tion information, afforded by measuring spatial vari- 
ations in the surface number density of faint back- 
ground galaxies (Broadhurst et al. 1995; Taylor et al. 
1998; Zhang & Pen 2005; Umetsu & Broadhurst 2008; 
Rozo & Schmidt 2010; Van Waerbeke et al. 2010), which 
we show here is readily detectable in the high-quality 
images we have obtained for this purpose. In our ear- 
lier work, particularly on A1689 and C10024+17, the 
magnification information determined from the back- 
ground counts was used only as a consistency check 
of distortion measurements (Broadhurst et al. 2008; 
Umetsu et al. 2010). Recently, we have shown in 
Umetsu & Broadhurst (2008) how to overcome the in- 
trinsic clustering of background galaxies, which other- 
wise perturbs locally the magnification signal, and how 
to combine the two independent weak-lensing data sets to 
improve the quality of two-dimensional mass reconstruc- 
tion using regularized maximum-likelihood techniques. 
Here we further explore new statistical methods designed 
to obtain an optimal combination of the complementary 
lens distortion and magnification effects. Our aim here is 
to develop and apply techniques to a sample of five well- 
studied, high-mass clusters with M > 10 15 M Q , A1689, 
A1703, A370, C10024+17, and RXJ1347-11, for examin- 
ing the underlying mass profiles extracted from the com- 
bined weak-lensing data sets. We then add to this de- 
tailed strong-lensing information for the inner £ 200 kpc 
region of these clusters, for which we have identified 
many new sets of multiple images from Advanced Cam- 
era for Surveys (ACS) observations (Broadhurst et al. 
2005b; Zitrin et al. 2009, 2010) with the Hubble Space 
Telescope (HST), to derive improved inner mass profiles 
for a full determination of the entire mass profiles of the 
five well-studied clusters. 

The paper is organized as follows. We briefly summa- 
rize in § 2 the basis of cluster weak gravitational lens- 
ing. In § 3 we present our comprehensive lensing method 
in a Bayesian framework for a direct reconstruction of 
the projected cluster mass profile from combined weak- 
lensing shape distortion and magnification bias measure- 
ments. In § 4 we apply our method to Subaru weak- 
lensing observations of five massive clusters to derive pro- 
jected mass profiles to beyond the cluster virial radius; 
we also combine our new weak-lensing mass profiles with 
inner strong-lensing based information from HST/ ACS 
observations to make a full determination of the entire 
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cluster mass profiles. Finally, summary and discussions 
are given in § 5. 

Throughout this paper, e use the AB magnitude sys- 
tem, and adopt a concordance ACDM cosmology with 
ro = 0.3, A = 0.7, and h = ^/(lOOkms^ 1 Mpc" 1 ) = 
0.7. Errors represent a confidence level of 68.3% (la) 
unless otherwise stated. 

2. BASIS OF CLUSTER WEAK LENSING 

The deformation of the image for a background source 
can be described by the Jacobian matrix A a p (ct, /3 = 
1,2) of the lens mapping. 9 The real, symmetric Jacobian 
A a p can be decomposed as A a p = (1— n)6 a p— T a p, where 
S a p is Kronccker's delta, k is the lensing convergence, and 
Y a p is the trace-free, symmetric shear matrix, 



r 



a/3 : 



+7i 72 
72 -7i 



(1) 



with 7 Q being the components of spin-2 complex gravita- 
tional shear 7 := 71+472- In the strict weak lensing limit 
where k, 7I -C 1, T a p induces a quadrupolc anisotropy of 
the background image, which can be observed from ellip- 
ticities of background galaxy images (Kaiser et al. 1995). 
The local area distortion due to gravitational lensing, or 
magnification, is given by the inverse Jacobian determi- 
nant, 

1 _ 1 
At ~ det~4 ~ (1 -k) 2 - | 7 | 2 ' ( ) 

where we assume subcritical lensing, i.e., det.4(0) > 0. 
The lens magnification \i can be measured from charac- 
teristic variations in the number density of background 
galaxies (Broadhurst et al. 1995, see also § 3.2). 

The lensing convergence K is a weighted projection of 
the matter density contrast along the line of sight (e.g., 
Jain et al. 2000) . For gravitational lensing in the cluster 
regime (e.g., Umetsu 2010), K is expressed as k(6) = 
E~; t E(0), namely the projected mass density E(0) in 
units of the critical surface mass density for gravitational 
lensing, defined as 



4nGD, 



max 



(J. 



D ds (z s ) 



D s (z s ) 



,(3) 



where D s , Dd, and Dd s are the proper angular diame- 
ter distances from the observer to the source, from the 
observer to the deflecting lens, and from the lens to the 
source, respectively, and (/3) = (Dd s /D s ) is the mean dis- 
tance ratio averaged over the population of source galax- 
ies in the cluster field. 

In general, the observable quantity for quadrupolc 
weak lensing is not the gravitational shear 7 but the com- 
plex reduced shear (see § 3.1), 



9(0) 



7(A) 
1 - k(G) 



(4) 



in the subcritical regime where det.4 > (or 1/g* in 
the negative parity region with det.4 < 0). The re- 

9 Throughout the paper we assume in our weak lensing analy- 
sis that the angular size of background galaxy images is sufficiently 
small compared to the scale over which the underlying lensing fields 
vary, so that the higher-order weak lensing effects, such as flex- 
ion, can be safely neglected; see, e.g., Goldberg & Bacon (2005); 
Okura et al. (2007, 2008). 



duced shear g is invariant under the following global lin- 
ear transformation: 



k(0) -> Xk(0) + 1 - A, 7(0) -> A 7 (0) 



(5) 



with an arbitrary scalar constant A 7^ 
(Schneider & Seitz 1995). This transformation is 
equivalent to scaling the Jacobian matrix .4(0) with 
A, A(6) — > XA(0), and hence leaves the critical curves 
det.4(0) = invariant. Furthermore, the curve k(0) = 1, 
on which the gravitational distortions disappear, is left 
invariant under the transformation (5). 

This mass-sheet degeneracy can be unambiguously 
broken by measuring the magnification effects (see § 3.2), 
because the magnification fi transforms under the invari- 
ance transformation (5) as 



M (0) -> A 2 M (0). 



(6) 



Alternatively, the constant A can be determined such 
that the mean k averaged over the outermost cluster re- 
gion vanishes, if a sufficiently wide sky coverage is avail- 
able. 10 

3. CLUSTER WEAK LENSING METHODOLOGY 

In this section we develop a Bayesian method to recon- 
struct the projected cluster mass profile E(0) from ob- 
servable lens distortion and magnification profiles, with- 
out assuming particular functional forms for the mass 
distribution, i.e., in a mo del- independent fashion. Al- 
though the methodology here is presented for the anal- 
ysis of individual clusters, it can be readily generalized 
for a statistical analysis using stacked lensing profiles of 
a sample of clusters. In § 4.6, alternatively, we provide a 
method to stack reconstructed projected mass profiles of 
individual clusters to obtain an ensemble-averaged pro- 
file. 

3.1. Lens Distortion Profile 

The observable quadrupolc distortion of an object due 
to gravitational lensing is described by the spin-2 re- 
duced shear, g = g\ + ig2 (equation [4]), which is co- 
ordinate dependent. For a given cluster center on the 
sky, one can form coordinate-independent quantities, the 
tangential distortion g + and the 45° rotated component, 
from linear combinations of the distortion coefficients as 
g + = — (g\ cos 2(f) + .92 sin 2(f)) and g x = — (172 cos 2(f) — 
gi sin 2<f>), with being the position angle of an object 
with respect to the cluster center. In the strict weak- 
lensing limit, the azimuthally-averaged tangential distor- 
tion profile g + w 7 + satisfies the following identity (e.g., 
Bartelmann & Schneider 2001): 7+(0) = «(< 6) - k{0), 
where k(9) is the azimuthal average of k(9) at radius 6, 
and k(< 9) is the mean convergence interior to radius 8. 
With the assumption of quasi-circular symmetry in the 
projected mass distribution (sec Umetsu & Broadhurst 
2008), the tangential distortion is expressed as 



9+(0) 



k(< 0) - k(6) 
1 - k(9) 



(7) 



10 Or, one may constrain the constant A such that the enclosed 
mass within a certain aperture is consistent with cluster mass es- 
timates from some other observations (e.g., Umetsu & Futamasc 
2000). 
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in the nonlinear but subcritical (det A(9) > 0) regime. 11 
In the absence of higher order effects, weak lensing only 
induces curl-free tangential distortions, while the az- 
imuthal averaged x component is expected to vanish. In 
practice, the presence of x modes can be used to check 
for systematic errors. 

From shape measurements of background galaxies, we 
calculate the weighted average of g + in a set of TV radial 
bands (i = 1, 2, TV) as 



g+m= E 



w 



E 



w (k) 9+(k) 



(8) 



where 9i is the center of the ith radial band of [0i, Oi+i], 
the index fc runs over all of the objects located within 
the ith annulus, <?+(&) is the tangential distortion of 
the fcth object, and wn^ is the statistical weight (see 
Umetsu k Broadhurst 2008; Umetsu et al. 2009, 2010) 
for the fcth object, given by 

w (k) = -3 ~, „. 2 ( 9 ) 



'g(k) 



at 



with a 2 ,, , being the variance for the shear estimate of the 
fcth galaxy and a 2 being the softening constant variance 
(e.g., Hamana et al. 2003). In our analysis, we choose 
a g = 0.4, which is a typical value of the mean rms o~ g 
over the background sample (e.g., Hamana et al. 2003; 
Umetsu & Broadhurst 2008; Umetsu et al. 2009, 2010). 
Wc use the continuous limit of the area-weighted band 
center for 9i (see equation [B3]). We perform a boot- 
strap error analysis to assess the uncertainty er+(#) in 
the tangential distortion profile g+(9). 

3.2. Magnification Bias Profile 

Lensing magnification, /i(0), influences the observed 
surface density of background sources, expanding the 
area of sky, and enhancing the observed flux of back- 
ground sources (Broadhurst et al. 1995). The former ef- 
fect reduces the effective observing area in the source 
plane, decreasing the number of background sources per 
solid angle; on the other hand, the latter effect amplifies 
the flux of background sources, increasing the number of 
sources above the limiting flux. The net effect is known 
as magnification bias, and depends on the intrinsic slope 
of the luminosity function of background sources. 

The number counts for a given magnitude cutoff m cu tj 
approximated locally as a power-law cut with slope s = 
dlog 10 TV(< m)/dm (s > 0), are modified in the presence 
of lensing as 



TV(< meat) « TV (< m cut )/-t 



2.5s-l 



(10) 



(Broadhurst et al. 1995), where TV (< m cut ) is the un- 
lensed counts, and \i is the magnification, /i = (1 — 
k)~ 2 (1 — \g\ 2 )~ 1 . In the strict weak lensing limit, the 
magnification bias is 5N/Nq ps (5s — 2)k. 

11 In general, a wide spread of the redshift distribution of back- 
ground galaxies, in conjunction with the single-plane approxima- 
tion, may lead to an overestimate of the gravitational shear in the 
nonlinear regime (Hoekstra et al. 2000). Thanks to the deep Sub- 
aru photometry, we found that this bias in the observed reduced 
shear is approximately Ag/g ~ 0.02ft to the first order of ft. See 
§ 3.4 of Umetsu et al. (2010) for details. 



For the number counts to measure magnification, wc 
use a sample of red background galaxies, for which 
the intrinsic count slope s at faint magnitudes is rel- 
atively flat, s ~ 0.1, so that a net count depiction 
results (Broadhurst et al. 2005a; Umetsu & Broadhurst 
2008; Umetsu et al. 2010). On the other hand, the faint 
blue background population tends to have a steeper in- 
trinsic count slope close to the lensing invariant slope 
(s = 0.4). The count-in-cell statistic TV(0; < m cut ) is 
measured from a flux-limited sample of red background 
galaxies on a regular grid of equal-area cells, each with 
a solid angle of Ail = (A9) 2 . Note that a practical 
difficulty of the magnification bias measurement is con- 
tamination due to the intrinsic clustering of background 
galaxies, which locally can be larger than the lensing- 
induced signal in a given cell. In order to obtain a clean 
measure of the lensing signal, such intrinsic clustering 
needs to be downweighted (e.g., Broadhurst et al. 1995; 
Zhang & Pen 2005). 

For a mass profile analysis, we calculate the mean num- 
ber density n fJi (9 i ) = dN(9i)/dil of the red background 
sample as a function of radius from the cluster center, 
by azimuthally averaging N{9), using the same radial 
bins (i = 1,2,..., TV) as done for the distortion mea- 
surement. The lens magnification bias is expressed in 
terms of the number density of background galaxies as 
77,^(0) = nofi(9) 2 - 5s ~ 1 with no being the unlcnsed mean 
surface number density of background galaxies. The nor- 
malization and slope parameters (no, s) can be estimated 
from the source counts in cluster outskirts using wide- 
field imaging data (§4.2). In practice, we adopt the fol- 
lowing prescription: 

• A positive tail of > 3tr cells is excluded in each 
annulus to remove inherent small scale clustering 
of the background (Broadhurst et al. 2008). 

• Each grid cell is weighted by the fraction of 
its area lying within the respective annular 
bins (Marshall et al. 2002; Umetsu & Broadhurst 
2008). 

• The uncertainty 0^(0) in 71^(9) includes not only 
the Poisson contribution but also the variance due 
to variations of the counts along the azimuthal di- 
rection, i.e., contributions from the intrinsic clus- 
tering of background galaxies (Umetsu et al. 2010). 

• The cell size A9 can be as large as the typical ra- 
dial band width for a mass profile analysis, which 
can cause an additional variance due to Poisson 
and sampling errors. We thus average over a set of 
radial profiles obtained using different girds offset 
with respect to each other by half a grid spacing in 
each direction. 

• The masking effect due to bright cluster galax- 
ies, bright foreground objects, and saturated pixels 
is properly taken into account and corrected for 
(Umetsu & Broadhurst 2008). In our analysis, we 
use Method B of Appendix A developed in this 
work. 

3.3. Bayesian Mass Profile Reconstruction 
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The relation between distortion and convergence is 
nonlocal, and the convergence derived from distortion 
data alone suffers from a mass-sheet degeneracy (§ 2). 
However, by combining the distortion and magnification 
measurements the convergence can be obtained unam- 
biguously with the correct mass normalization. Here we 
aim to derive a discrete convergence profile from observ- 
able lens distortion and magnification profiles (see § 3.1 
and § 3.2) within a Bayesian statistical framework, al- 
lowing for a full parameter-space extraction of model 
and calibration parameters. A proper Bayesian statis- 
tical analysis is of particular importance to explore the 
entire parameter space and investigate the parameter de- 
generacies, arising in part from the mass-sheet degener- 
acy. 

In this framework, we sample from the posterior prob- 
ability density function (PDF) of the underlying signal 
s given the data d, P(s\d). Expectation values of any 
statistic of the signal s shall converge to the expectation 
values of the a posteriori marginalized PDF, P(s\d). In 
our problem, the signal s is a vector containing the dis- 
crete convergence profile, k, = n{0i) (i = 1,2, ..,N), and 
the average convergence within the inner radial boundary 
6 m i a = Q\ of the weak lensing data, 7c m i n = k(< m in), so 
that s = {K m i n , fti, «2, «at}, being specified by (N + 1) 
parameters. The Bayes' theorem states that 

P{s\d) oc P(s)P(d\s), (11) 

where £(s) = P(d\s) is the likelihood of the data given 
the model (s), and P(s) is the prior probability distri- 
bution for the model parameters. 

3.3.1. Weak Lensing Likelihood Function 

We combine complementary and independent weak 
lensing information of tangential distortion and magnifi- 
cation bias to constrain the underlying cluster mass dis- 
tribution s = {/Cmin, Kj}i=i- The total likelihood func- 
tion C of the model s for combined weak lensing obser- 
vations is given as a product of the two separate likeli- 
hoods, C = C g Cfj_, where C g and £ M are the likelihood 
functions for distortion and magnification, respectively. 
The log-likelihood for the tangential distortion is given 
as 

- ln£ 3 = \ J2 f+»] 2 + const., (12) 

2 i=i a +-* 

where g+ t i(s) is the theoretical prediction for the ob- 
served distortion g+j, and the errors (i = 1, 2, N) 
due primarily to the variance of the intrinsic source cllip- 
ticity distribution can be conservatively estimated from 
the data using Bootstrap techniques (§ 3.1). Similarly, 
the log-likelihood function for the magnification bias is 
given as 

1 x—^ \n„ i — h.j As)} 2 

- In£ M = - 2 + const -' ( 13 ) 

i=i °W 

where h^^(s) is the theoretical prediction for the ob- 
served counts n^j, and the errors a^.i include both con- 
tributions from Poisson errors in the counts, 0p o i ssorl) j, 
and contamination due to intrinsic clustering of red back- 



Table 1 

The Cluster Sample: Rcdshift and Subaru Filter Information 



Cluster Redshift Filters Detection band Seeing a Refs. 

(arcsec) 



A1689 


0.183 


i'V 


i' 


0.82 


1 


A1703 


0.281 


g'r'i' 


r' 


0.78 


2.3 


A370 


0.375 




Rc 


0.60 


2.3 


C10024+17 


0.395 


BjR c z' 


Rc 


0.80 


4 


RXJ1347-11 


0.451 


VR c z' 


Rc 


0.76 


2.3 



References. [1] Umetsu & Broadhurst (2008); [2] 

Broadhurst et al. (2008); [3] Medezinski et al. (2010); [4] Umetsu et al. 
(2010). 

Note. — For observational details, see the references. Refs. [2] and 
[3] adopted a cluster redshift of z = 0.258 for A1703. 
a Seeing FWHM in units of arcsec in the final co-added detection image. 



Table 2 

Einstein Radius Information 



Cluster 


Einstein radius, 9 e j n 


Refs. 




(arcsec) 




A1689 


53 ± 3" (z s = 3.04) a 


1,2 


A1703 


31 ± 3" (z s = 2.627) 


3, 4, 5 


A370 


37 ± 3" (z s = 2) 


6, 7 


C10024+17 


30 ±3" (z s = 1.675) 


8 


RXJ1347-11 


35 ± 2" (z s = 2.2) 


7, 9, 10 



References. — [1] Broadhurst et al. 
(2005b); [2] Coe et al. (2010); [3] 
Limousin et al. (2008); [4] Richard et al. 

(2009) ; [5] Zitrin et al. (2010); [6] Richard et al. 

(2010) ; [7] Zitrin et al. (2011), in preparation; 
[8] Zitrin et al. (2009); [9] Halkola et al. (2008); 
[10] Bradac et al. (2008). 

a This Einstein radius constraint translates into 
6» ein ~ 47" at z a = 2 (cf. 6» cin = 47.0 ± 1.2" at 
z s = 2 by Ref. [2]). 

ground galaxies, <7 c fust,i) as discussed in § 3.2: 

a H,i — CT Poisson,i + ^clust.r 

The lensing observables and n^; (i = 1,2,..., N) can 
be readily expressed by the given model parameters s, as 
shown in Appendices B and C. 

3.3.2. Prior Lnformation 

For each parameter of the model s, we consider a sim- 
ple flat prior with a lower bound of s = 0, that is, 

Kmin>0, (15) 

k 4 >0 (i = l,2,...,A0- (16) 

Additionally, we account for the calibration uncertainty 
in the observational parameters, i.e., the normalization 
and slope parameters (no,s) of the background counts 
and the relative lensing depth uj due to population-to- 
population variations between the background samples 
used for the magnification and distortion measurements 
(see Appendix C). 

4. APPLICATIONS: SUBARU OBSERVATIONS OF FIVE 
STRONG-LENSING CLUSTERS 

4.1. Cluster Sample and Observations 
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Table 3 

Background galaxy samples 



Cluster 


Distortion analysis (full background) 
(arcmin 2 ) 


Magnification analysis (red background) 

™ z a ,g </9> 
(arcmin -2 ) 




A1689 a 


8.8 


0.71 ±0.12 


0.69 


12.0 


0.71 ±0.12 


0.69 


1.00 


A1703 


10.0 


1.10 ±0.18 


0.68 


6.9 


0.93 ±0.12 


0.64 


0.95 


A370 b 


16.7 


1.11 ±0.12 


0.58 


21.6 


1.05 ±0.11 


0.57 


0.98 


C10024+17 


17.2 


1.29 ±0.15 


0.61 


17.8 


1.14 ±0.09 


0.56 


0.92 


RXJ1347-11 


6.4 


0.98 ± 0.06 


0.47 


7.7 


1.04 ±0.07 


0.49 


1.05 



Note. — Wc assume 5% uncertainty in our estimates of the relative lensing depth uj. 
a Mean surface number density of background galaxies. 

Effective mean source redshift of the background sample corresponding to the mean depth (/3), defined such 
that &(z 8 ,&) = 03). 

c Distance ratio averaged over the redshift distribution of the background sample, taken from 
Umetsu & Broadhurst (2008, A1689), Medezinski et al. (2010, A1703, A370, RXJ1347-11) and Umctsu ct al. 
(2010, C10024+17). 

Relative lensing depth of the red background sample, uj — (/3(red))/(/?(full)), with respect to the full 
background sample. 



Table 4 

Normalization and slope of red background counts 



Cluster 


m cut a 


Normalization, no 


Slope, s c 




(AB mag) 


(arcmin -2 ) 




A1689 


i' = 25.5 


11.5 ±0.3 


0.124 ±0.060 


A1703 


i' = 26 


8.2 ±0.2 


0.099 ± 0.070 


A370 


z' = 26 


22.5 ±0.4 


0.042 ± 0.070 


C10024+17 


z' = 25.5 


18.3 ±0.3 


0.121 ± 0.040 


RXJ1347-11 


z' = 26 


7.0 ±0.3 


0.025 ± 0.070 


a Fainter magnitude cut of the red background sam 


pie, m cut . 


k Normalization of unlensed red background counts 


no. 


c Slope of the 


unlensed red 


background counts 


at m cu t , s — 


dlog 10 N (< m c 


at ) /dm. 







In this section wc apply our Bayesian mass reconstruc- 
tion method to a sample of five well-studied high-mass 
clusters (M > 10 15 M Q ) at intermediate redshifts, A1689 
[z = 0.183), A1703 {z = 0.281), A370 (z = 0.375), 
C10024+17 (z = 0.395), and RXJ1347-11 (z = 0.451), 
observed with the wide- field camera (34' x 27'; see 
Miyazaki et al. 2002) on the 8.2m Subaru telescope. Ta- 
ble 1 gives a summary of cluster observations. The clus- 
ters were observed deeply in several optical passbands, 
with exposures in the range 2000-10000 s per passband, 
with seeing in the detection images ranging from 0.6" to 
0.8" (see Umetsu & Broadhurst 2008; Broadhurst et al. 
2008; Medezinski et al. 2010; Umetsu et al. 2010). 

These massive clusters are known as strong lens- 
ing clusters, displaying prominent strong-lensing fea- 
tures and large Einstein radii of 6* e ; n ^ 30" (e.g., for a 
fiducial source redshift z s ~ 2; Broadhurst & Barkana 
2008; Oguri & Blandford 2009). For the clusters, 
the central mass distributions have been recovered in 
detail by our strong-lens modeling (Broadhurst et al. 
2005b; Zitrin et al. 2009, 2010). Here the models 
were constrained by a number of multiply-lensed im- 
ages identified previously in very deep multi-color 
imaging with HST/ ACS (e.g., Broadhurst et al. 2005b; 
Bradac et al. 2008; Halkola ct al. 2008; Limousin et al. 
2008; Richard et al. 2009, 2010; Zitrin et al. 2009, 2010). 
Table 2 gives a summary of the Einstein radii of the 
five clusters as constrained from our detailed strong lens 



modeling with the ACS observations. 

4.2. Background Galaxy Samples 

Here wc will add to our high-quality spin-2 
shape measurements (Umetsu & Broadhurst 2008; 
Medezinski et al. 2010; Umetsu et al. 2010) the in- 
dependent magnification information based on deep 
multi-band imaging with Subaru, in order to achieve the 
maximum possible lensing precision. Full details of the 
Subaru observations and weak-lensing shape analysis 
of these clusters were presented in a series of our 
papers (Umetsu & Broadhurst 2008; Broadhurst et al. 
2008; Umetsu et al. 2009; Medezinski et al. 2010; 
Umetsu et al. 2010; Zitrin et al. 2010, see also Table 
1). The level of shear calibration bias with our imple- 
mentation of the KSB+ method (Kaiser et al. 1995) has 
been assessed by Umetsu et al. (2010) using simulated 
Subaru Suprime-Cam images (M. Oguri 2010, in private 
communication; Massey et al. 2007). We find, typically, 
|m| ^ 5% of the shear calibration bias, and c ~ 10~ 3 
of the residual shear offset which is about 1 order of 
magnitude smaller than the typical distortion signal 
in cluster outskirts. This level of calibration bias is 
subdominant compared to the statistical uncertainty 
(AM/M ~ 15%) due to the intrinsic scatter in galaxy 
shapes, and to the dilution effect which can lead to an 
underestimation of the true signal for R ^ 400kpc/i~ 1 
by a factor of 2-5 (see Figure 1 of Broadhurst et al. 
2005a). 

A careful background selection is critical for a weak- 
lensing analysis so that unlensed cluster members 
and foreground galaxies do not dilute the true lens- 
ing signal of the background (Broadhurst et al. 2005a; 
Medezinski et al. 2007, 2010). We use undiluted samples 
of background galaxies derived in our previous lensing 
work, as summarized in Table 3. When deep multi- 
color photometry is available in our cluster fields, we 
use the background selection method of Medezinski et al. 
(2010, see also Umetsu ct al. 2010) to define blue 
and red background samples (A1703, A370, C10024+17, 
RXJ1347-11), which relies on empirical correlations for 
galaxies in color-color space derived from the deep Sub- 
aru photometry, by reference to the deep photometric- 
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Figure 1. Radial profiles of the averaged spin-2 distortion and magnification bias measurements based on background galaxies registered 
in deep Subaru images, shown separately for the five clusters, A1689 (top-left), A1703 (top-middle), A370 (top-right), C10024+17 (bottom- 
left), and RXJ1347-11 (bottom-right). For each cluster, the top panels show the radial profiles (open squares and error bars) of the 
tangential distortion (upper panel) and the 45° rotated (x) component (lower panel) based on the spin-2 shape measurements of the full 
background galaxy sample (Table 3). The bottom panels show the count profiles of red background galaxies. The squares and triangles 
(upper panel) show the respective results with and without the mask correction due to bright foreground objects and cluster members. The 
blue hatched area represents the 68.3% confidence bounds for the predicted count depletion curve from an NFW model constrained by our 
Subaru distortion analysis, demonstrating clear consistency between these two independent lensing obscrvablcs. The gray horizontal bar 
represents the constraints on the unlensed count normalization no. The histogram in the lower panel shows the observed counts of the red 
background galaxies in each annular bin. 



rcdshift survey in the COSMOS field (Ilbcrt et al. 
2009). Otherwise, we use the color-magnitude selection 
method (Broadhurst et al. 2005a; Medezinski et al. 2007; 
Umetsu & Broadhurst 2008) to define a sample of red 
galaxies (A1689) whose colors are redder than the red 
sequence of cluster E/SO galaxies. These red galaxies 
are expected to lie in the background by virtue of K- 
corrections which are greater than the red cluster se- 
quence galaxies, as convincingly demonstrated spectro- 
scopically by Rines & Geller (2008). Apparent magni- 
tude cuts are applied in the reddest band available for 
each cluster to avoid incompleteness near the detection 
limit. 



A flux-limited sample of red background galaxies is 
used for the magnification analysis (see § 3.2). For 
the distortion analysis, we use a full composite sam- 
ple of red and blue (if available) background galaxies, 
where the galaxies used are well resolved to make reli- 
able shape measurements (see, e.g., Umetsu et al. 2010). 
Table 3 lists for respective color samples the mean sur- 
face number density n of background galaxies, the effec- 
tive source redshift, z s „a, 12 and the mean distance ra- 

12 The effective single-lens plane redshift z s ,Di corresponding to 
the mean depth (/3), is defined as (0) = /3fz s a). For details, see 
§ 3.4 of Umetsu et al. (2010). 



8 



Bayesian Analysis of Weak Lensing Distortion and Magnification Measurements 




Figure 2. One-dimensional marginalized posterior probability density functions (PDFs) for the discrete mass profile s = 
(Ttjnin, «i, K2, Kn), shown for A1689 (see Figure 3). The solid vertical lines show the peak locations (mode values) of the a poste- 
riori marginalized distribution for each parameter. The results arc marginalized over all other parameters, including the observational 
parameters (no,s,uj). The resulting posterior distributions are all single-peaked, and nearly Gaussian for most of the model parameters. 
The mass-sheet degeneracy is broken thanks to the inclusion of magnification data. 
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Figure 3. Reconstructed projected mass profiles k(8) = £(c9)/E cr it for the five clusters, A1689 (top left), A1703 (top middle), A370 
(top right), C10024+17 (bottom left), and RXJ1347-11 (bottom right). The open boxes show the results obtained using the Bayesian 
mass reconstruction method from the combined Subaru tangential distortion and magnification bias measurements of background galaxies 
(Figure 1). The innermost (first) box represents the average convergence 7t(< mlrl ) interior to the minimum radius 8 lrlln of the weak lensing 
data. Also shown for each cluster is the best-fit NFW profile (solid curve) for the combined Subaru distortion and magnification data. The 
gray boxes represent the mass profile reconstructed using the nonlinear extension of aperture mass densitometry by Umetsu & Broadhurst 
(2008) based on the same tangential distortion data (but without the magnification data combined), which employs an outer boundary 
condition on the mean convergence in the outermost radial bin. The horizontal extent of each box represents the range of the radial band, 
and its vertical extent shows the statistical 1<t uncertainty in ft. The cross in each open box indicates the band-averaged convergence, being 
marked at the area-weighted center of the radial band. The errors are correlated in both reconstructions. 



tio (/3) averaged over the source redshift distribution. 
We also quote in Table 3 the values of the relative 
mean depth, oj = (/3(red))/(/3(full)}, between the back- 
ground samples used for the magnification and distortion 
measurements (see § 3.3.2 and Appendix C). Wc adopt 
a conservative uncertainty of 5% in the relative mean 
depth u> for all of the clusters based on our previous 
work (Umetsu k Broadhurst 2008; Umetsu et al. 2009; 
Medezinski et al. 2010; Umetsu et al. 2010). 

The conversion from the observed counts into magni- 
fication depends on the normalization no and the slope 
parameter s of the unlensed number counts, which can 
be reliably estimated thanks to the wide field of view 



of Subaru/Suprime-Cam (see § 4.4). Table 4 lists the 
magnitude cuts (m cut ) and unlensed count parameters 
(no,s) as measured from the red background counts in 
outermost annular regions in cluster outskirts ( ^ 10'). 

4.3. Subaru Weak Lensing Profiles 

Following the methodology outlined in § 3.1 and § 3.2 
(see also Appendix A), we derive lens distortion and mag- 
nification profiles of five massive clusters from Subaru ob- 
servations. In order to obtain meaningful radial profiles, 
one must carefully define the center of the cluster. It is 
often assumed that the cluster mass centroid coincides 
with the position of the brightest cluster galaxy (BCG), 



Umetsu et al. 
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Table 5 

Subaru weak lensing data 



Cluster 


^min 3 Pmax 


N 




Detection S/N 




(') 




9+ 


rip. 


Combined (ft) 


A1689 


1, 18 


11 


13.8 


8.8 


17.8 


A 1703 


1, 16 


10 


9.9 


7.1 


12.7 


A370 


0.7, 16 


14 


17.3 


11.0 


23.8 


C10024+17 


0.6,15 


12 


13.0 


8.5 


18.5 


RXJ1347-11 


0.8, 11 


11 


9.7 


5.7 


12.6 



Note. — The lensing profiles are calculated in N discrete radial 
bins over the radial range of 9 — [$ m i n , $ max ], with a logarithmic 
radial spacing of A ln# — ln(# max /# m i n )/iV. For each cluster, # m i n 
is taken such that # m i n > 6 e \ n (see Table 2). 



whereas the BCGs can be offset from the mass centroids 
of the corresponding dark matter halos (Johnston et al. 
2007; Oguri et al. 2010; Oguri & Takada 2011). Here we 
utilize our detailed mass maps in the cluster cores re- 
covered from strong-lens modeling of ACS observations 
(§4.1), providing an independent mass-centroid determi- 
nation. We find that for these five clusters there is only 
a small offset of typically £j 5" (20kpc/i -1 at the high- 
est cluster redshift of our study, Zd = 0.451) between the 
BCG and the dark-matter center of mass (see also § 4.2 
of Umetsu et al. 2010), often implied by other massive 
bright galaxies in the vicinity of the BCG. This level of 
cluster centering offset is substantially small as compared 
to the typical inner radial boundary of weak lensing ob- 
servations, f? m m ~ 1'. In the present work, we therefore 
simply assume that the cluster mass centroid coincides 
with the location of the BCG, which is adopted as the 
cluster center in our one-dimensional profile analysis. 

The lensing profiles were calculated in N discrete 
radial bins over the range of radii 9 = [9 m - m , f? max ] , 
with a constant logarithmic radial spacing of A In 9 = 
ln(9 max /6 m - m )/N, where the inner radial boundary 6? m i n 
is taken such that 9 m i n > 9 cm (sec Table 2). The typi- 
cal inner boundary is 9 mm ~ l'(> f? c i n ) for cluster weak 
lensing. The outer radial boundary f? max was chosen to 
be sufficiently larger than the typical virial radius of high 
mass clusters with M vir = (1 - 2) x 10 15 M Q h -1 , i? max = 
Ai#max ~ 3 Mpc, but sufficiently small (9 max £ 18') with 
respect to the size of the Suprime-Cam's field-of-view 
so as to ensure accurate PSF (point spread function) 
anisotropy correction. The number of radial bins N was 
determined for each cluster such that the pcr-pixel detec- 
tion signal-to-noisc ratio (S/N) is of the order of unity, 13 
which is optimal for an inversion problem. The radial 
binning scheme is summarized in Table 5. 

In Figure 1 we show the resulting distortion and mag- 
nification profiles for our sample of five massive lensing 
clusters. In all the clusters, a strong depletion of the 
red galaxy counts is seen in the central, high-density re- 
gion of the cluster, and clearly detected out to several 
arcminutes from the cluster center. The statistical sig- 
nificance of the detection of the depletion signal is in the 
range 6cr-ll<7 (see Table 5). The detection significance 
of the tangential distortion derived from the full back- 
ground sample ranges from 10cr to 17cr, and is better 
than the red counts (see, e.g., Bartclmann & Schneider 
1999). The x -component is consistent with a null signal 

13 We quantify the significance of a detection for a given lensing 
profile in analogy to equation (38) of Umetsu & Broadhurst (2008). 



in most of radial bins, indicating the reliability of our 
wcak-lcnsing analysis. The magnification measurements 
with and without the masking correction are roughly con- 
sistent with each other. Typically, the masking area is 
negligible (a few %) at large radii, and increases up to 
(10 — 20)% of the sky close to the cluster center (sec 
Appendix A). To test the consistency between our dis- 
tortion and depletion measurements, we calculate the de- 
pletion of the counts expected for the best-fitting NFW 
profile derived from the distortion measurements (Fig- 
ure 1), normalized to the observed density no (Table 4). 
This comparison shows clear consistency between two in- 
dependent lensing obscrvablcs with different systcmatics 
(see § 5.5 of Umetsu & Broadhurst 2008), which strongly 
supports the reliability and accuracy of our weak-lcnsing 
analysis (see also supplemental material presented in Ap- 
pendix D). 

4.4. Cluster Mass Profile Reconstruction 

We use a Markov chain Monte Carlo (MCMC) ap- 
proach with Metropolis-Hastings sampling to reconstruct 
the discrete cluster mass profile s = {7« mm , Ki}fL 1 within 
a Bayesian statistical framework (§ 3.3). We largely fol- 
low the sampling procedure outlined in Dunklcy et al. 
(2005), but employ the Gclman- Rubin R statistic 
(Gelman & Rubin 1992) as a simple but reasonable con- 
vergence criterion of generated chains. Once convergence 
to a stationary distribution is achieved, we run a long 
final chain of 300,000 steps, which adequately samples 
the underlying posterior probability distribution. For all 
of the parameters, the number of iterations required for 
convergence is much less than our final chain length. 
Note, only the final chain is used for our parameter 
estimations and error analysis. We estimate the loca- 
tion of the maximum a posteriori probability for each 
model parameter using the bisection method in conjunc- 
tion with bootstrap techniques. The covariance matrix 
Cij = 1, 2, N + 1) for the discrete mass profile s 
is estimated from the MCMC samples. As an example, 
we show in Figure 2 one-dimensional marginalized pos- 
terior PDFs for the mass profile s of A1689 (N = 11). 
The results are marginalized over all other parameters, 
including the observational parameters (no,s,u>). 

The resulting posterior distributions are all clearly 
single-peaked, and approximately Gaussian for most of 
the parameters. It is clearly evident that the mass-sheet 
degeneracy is broken thanks to the inclusion of magnifi- 
cation information on the local area distortion. Exclud- 
ing magnification data, on the other hand, strongly mod- 
ifies the Gaussian shape of the marginalized posterior 
PDFs, producing long non-Gaussian tails and broaden- 
ing the distribution function, resulting in large errors for 
the reconstructed mass profile. The improvement here 
from adding the magnification measurements is signifi- 
cant, ~ 30% in terms of cluster mass profile measure- 
ments (see Table 5). 

Figure 3 shows the lensing convergence profiles s = 
{^min, ^1, k,2, k-n}, for the five clusters reconstructed 
using our Bayesian method from combined Subaru dis- 
tortion and magnification data. Also shown for com- 
parison are independent reconstructions from the same 
tangential distortion data (but without the magnifica- 
tion data combined) using the one-dimensional method 
of Umetsu & Broadhurst (2008, see also Umetsu et al. 
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Table 6 

Best-fit NFW model parameters 



Cluster 


M viT 
(lO^Moh" 1 ) 


NFW (weak lensing) 

c v ir X 2 /dof 


n a 

(") 


M vlr 
(10 15 M Q h- 1 ) 


gNFW (weak+strong lensing) 
c_2 b ol x 2 /dof 


n a 
y cin 

(") 


A1689 
A1703 
A370 
C10024+17 
RXJ1347-11 


-, 9S9 +0.217 
l.zoz_ 176 

-, r,or, + 0.244 

l.zoz_ Q 2()4 

2 451+0-309 
z -<± d -0.262 

1 o , 7c"r"0.232 
±.oto_ 2Q1 

1 48S+0- 277 
1.4B8_ 244 


12.80+J^ 

7 rio + 2.36 
'• uz -1.70 

7 00+1.09 

'• uu -0.92 
»-»2_ 1 gg 

9.08+j; 1 ^ 


4.2/10 
6.1/9 
10.2/13 
11.7/11 
10.4/10 


48.7lf 3 ;l 

28-i±ii:i 

50 4+ 11 - 
36 q+ 113 

■3D-»_10.1 

43 2+ 13 6 
40 -12.6 


1 301+ 0193 

i 979 +0.234 
-0.191 
9 9 y^+0-255 
z.zio_ 220 
1 QQO+0-247 
l-JJ»_0.203 

1 435+0.H4 


7 07+1-08 
'• u '-1.06 

8-05±i-« 

7 16+ ' 43 
'• 1D -0.43 


0.268t°;«| 
0.934«;^ 

0.392+S 
0.814+°-??? 
046+° :205 


4.4/20 
7.2/21 
15.2/26 
11.1/23 
58.2/28 


49.7+^ 2 2 9 

27 9+ 14 - 7 
z '- a -14.0 

30 5+ 9 ' 4 

oi 0+13.4 

0± -°-9.2 


a Einstein radius 
K 


arcdictcd by the be 


st-fit NFW (gNFW) model, ev 


aluatcd at the 


arc redshift given in 


Table 2. 









° Effective concentration parameter for gNFW, c_ 2 = ^ v i r / r -2 = c vir/( 2 ~ Q )- 

2010) based on the nonlinear extension of aperture mass 
densitometry, which employs an outer boundary con- 
dition on the mean convergence in the outermost ra- 
dial bin, ?? max = K(0jv-i,0max)- Here K max for an iso- 
lated NFW halo can be negligibly small if the outer- 
most radii are taken as large as the cluster virial ra- 
dius (see Umetsu et al. 2010). This method has been ap- 
plied successfully to Subaru weak lensing observations of 
massive clusters including A1689 (Umetsu & Broadhurst 
2008; Umetsu et al. 2009), A1703 (Oguri et al. 2009; 
Zitrin et al. 2010), and C10024+17 (Umetsu et al. 2010). 
Our results with different combinations of lensing mea- 
surements and boundary conditions, having different sys- 
tematics, are in agreement with each other. This con- 
sistency clearly demonstrates that our results are ro- 
bust and insensitive to the choice of boundary condition 
as well as to systematic errors in the lensing measure- 
ments, such as the shear calibration error, as found by 
Umetsu & Broadhurst (2008). 

Unlike the distortion effect, the magnification bias 
due to the local area distortion falls off sharply with 
increasing distance from the cluster center. We find 
from the reconstructed mass profiles that the lens con- 
vergence at large radii of 9 = [10', 15'] is of the order 
K = 5 X 10~ 3 - 0.01. The expected level of the deple- 
tion signal in the weak-lensing limit is Sh^/hq m —2k 
for a mtmmaZZy-dcpleted sample with s = 0, indicat- 
ing a depletion signal of ;$ (1 — 2)% in the cluster out- 
skirts where we have estimated the unlensed background 
counts, uq. This level of signal is smaller than the frac- 
tional uncertainties in estimated unlensed counts no of 
(2 — 4)%, thus consistent with the assumption. Note that 
the calibration uncertainties in our observational param- 
eters (no, s, u>) have been properly taken into account and 
marginalized over in our Bayesian analysis. 

In the presence of magnification, one probes the num- 
ber counts at an effectively fainter limiting magnitude: 
?ricut + 2.5 log 10 u(9). The level of magnification is on 
average small in the weak-lensing regime but for the in- 
nermost bin reaches a factor of 2 to 4 depending on the 
cluster. Here we use the count slope at the fainter ef- 
fective limit (mum) when making the magnification es- 
timate, to be self-consistent. In our analysis we have 
implicitly assumed that the power-law behavior (equa- 
tion [10]) persists down to ~ 1 mag fainter than m cu t 
where the count slope may be shallower. For a given 
level of count depletion, an ovcrestimation of the count 
slope could lead to an overestimation of the magnifica- 
tion, thus biasing the resulting mass profile. However, 



the number count slope for our data flattens only slowly 
with depth varying from s ~ 0.1 to s ~ 0.05 from a limit 
of m = 25.5 to to = 26.5, so that this introduces a small 
correction of only typically 8%— 11% for the most mag- 
nified bins (u = 2 — 4). In fact, we have found a good 
consistency between the purely shear-based results and 
the results based on the combined distortion and magni- 
fication data (see Figure 3). 

To quantify and characterize the cluster mass distri- 
bution, we compare the reconstructed k profile with 
the physically and observationally motivated NFW 
model. Here we consider a generalized paramctrization 
of the NFW model of the following form (Zhao 1996; 
Jing & Suto 2000): 



p{r) 



(r/r s ) a (l + r/r s 



(17) 



which has an arbitrary power-law shaped central cusp, 
73D = —ot, and an asymptotic outer slope of 73D = —3. 
This reduces to the NFW model for a = 1. We refer 
to the profile given by equation (17) as the generalized 
NFW (gNFW, hereafter) profile. It is useful to introduce 
the radius r_2 at which the logarithmic slope of the den- 
sity is isothermal, i.e., 73D = —2. For the gNFW profile, 
r_2 = (2 — a)r s , and thus the corresponding concentra- 
tion parameter reduces to c_2 = r v i r /r_2 = c V i r /(2 — a). 
We specify the gNFW model with the central cusp slope, 
a, the halo virial mass, M vir , and the concentration, 
c_2 = c v i r /(2 — a). We employ the radial dependence 
of the gNFW lensing profiles given by Keeton (2001). 

We first fix the central cusp slope to a = 1 (NFW), 
and constrain (M v i r ,c v i r ) from % 2 fitting to the discrete 
cluster mass profile s = {7t m hi, «i}£Li reconstructed from 
the combined weak-lensing distortion and magnification 
measurements. The y 2 function for weak lensing is de- 
fined by 



JV+l JV+1 

EE 



i(M 



vir j L-vir 



f (M, 



vm <--vir 



)]: 



(18) 

where s(M vu -, c V i r ) is the NFW model prediction for the 
discrete mass profile s. The resulting constraints on the 
NFW model parameters and the predicted Einstein ra- 
dius 9 c i D arc shown in Table 6. For all the cases, the 
best-fit NFW model from weak lensing properly repro- 
duces the observed location of the Einstein radius, con- 
sistent with the independent strong-lensing observations 
(see Table 2). In Table 5 we quote the values of the total 
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detection S/N in the reconstructed mass profile s based 
on the combined distortion and magnification data. 

4.5. Cluster Mass Estimates 

Our comprehensive Bayesian analysis of the weak lens- 
ing distortion and magnification of background galaxies 
allows us to recover the mass normalization, given as 
the mean convergence K m j n within the innermost mea- 
surement radius 9 mnl (> em ), without employing inner 
strong lensing information. 

We use the non-parametric deprojection method of 
Broadhurst & Barkana (2008) to derive for each cluster 
three-dimensional virial quantities (r v i r ,M v ; r ) and val- 
ues of Ma = M3d(< ?*a) within a sphere of a fixed 
mean interior overdensity A with respect to the critical 
density p cr it(zd) of the universe at the cluster redshift 
Zd- We first deproject the two-dimensional mass pro- 
files obtained in § 4.4 and derive non-parametric three- 
dimensional mass profiles M^(< r) simply assuming 
spherical symmetry, following the method introduced by 
Broadhurst & Barkana (2008). This method is based on 
the fact that the surface-mass density S(i?) is related to 
the three-dimensional mass density p(r) by an Abel inte- 
gral transform; or equivalently, one finds that the three- 
dimensional mass M$d{< t) out to spherical radius r is 
written in terms of £(i?) as 

M 3D {<r) = M 2D (< R min ) (19) 
+ 2tt J dRKE(R)-4:j dRRf (J^J E(iJ), 

where f(x) = (x 2 - l)" 1 / 2 - tan" 1 ^ 2 - l)" 1 / 2 
(Broadhurst & Barkana 2008; Umetsu et al. 2010), 14 
and the first term of the right-hand side can be obtained 
as M 2D (< -Rmin) = 7ri? mm £ cri t7cmin with R min = D d 9 min . 
The errors are estimated from Monte Carlo simulations 
based on the full covariance matrix of the lensing con- 
vergence profile (for details, see Umetsu et al. 2010). In 
Figure 4 we show the non-parametric cumulative pro- 
jected (blue-hatched) and spherical (gray-shaded) mass 
profiles, M 2 d(< 0) and M^u(< 0), separately for the five 
clusters. 

Table 7 gives a summary of the spherical mass es- 
timates M3d(< r&) corresponding to A = A v i r ,200, 
and 2500 from our non-parametric deprojection analy- 
sis, where A v i r ~ 110 — 130 is the virial overdensity 
of the spherical collapse model evaluated at the clus- 
ter redshift za- Overall, we find a good agreement be- 
tween the virial mass estimates from the parametric and 
non-parametric deprojection approaches (see § 4.7 for the 
results of gNFW fits to the combined weak and strong 
lensing data). 

Our virial mass estimate of A1689 is obtained as 
M vir = (1.30 ± 0.21) x lO^Mg/i' 1 from the com- 
bined distortion and magnification profiles, consistent 
with the results of Umetsu & Broadhurst (2008), who 
combined strong lensing, weak lensing distortion and 
magnification data in a full two-dimensional analysis, 
and derived M vir = l.5ta.t x 10 1B M©/i _1 , where this 

14 This integral transformation has an intcgrablc singularity 
at the lower limit of the second integrand (R = r), which can be 
readily avoided by a suitable coordinate transformation. 



la error includes both statistical and systematic uncer- 
tainties (see also Umctsu et al. 2009; Kawaharada et al. 
2010; Molnar et al. 2010). 15 This is also in good agree- 
ment with the results from the recent high-resolution 
lensing observations by Coe et al. (2010), M V { T = 
(lAtoz) x 10 15 Af o /i _1 . These lensing results are con- 
sistent with careful dynamical work by Lemze et al. 
(2009a), who obtained a virial mass estimate of M w - lr = 
(1.3 ± 0.4) x 1O 15 M h' 1 for A1689. As found early 
by Broadhurst et al. (2008), our comprehensive weak- 
lensing analysis implies A370 is the most massive cluster 
now known, M vir = (2.40 ± 0.25) x 1O 15 M /i _1 . Our 
virial mass estimate is slightly higher than that derived 
in our earlier weak-lensing work combined with the in- 
ner Einstein radius information, Af v ; r = (2.1 ± 0.2) x 
IO^Mq/i" 1 (Broadhurst et al. 2008), where the differ- 
ence is primarily due to our improved background selec- 
tion and depth estimate presented in Medezinski et al. 
(2010). For RXJ1347-11, we find its NFW virial mass is 
slightly overestimated compared to the non-parametric 
estimate of M vir = (1.15±0.25) x 10 15 M© bT 1 due to the 
projection of subclumps associated with the large scale 
structure around the cluster (see Lu et al. 2010). Our 
non-parametric mass estimates are in good agreement 
with independent X-ray, dynamical and lensing analyses 
(Kling et al. 2005; Miranda et al. 2008; Lu et al. 2010). 
For A1703 and C10024+17, our new mass estimates are 
fully consistent with our recent weak-lensing results de- 
rived from the Subaru distortion data alone (Zitrin et al. 
2010; Umctsu ct al. 2010, see also § 4.4). 

4.6. Stacking Analysis 

The statistical precision of lensing constraints can be 
further improved by stacking the signal from an ensemble 
of clusters with respect to their centers, providing aver- 
age properties of cluster mass profiles. As discussed by 
Okabe et al. (2010), this stacking analysis has several im- 
portant advantages. A notable advantage of the stacking 
analysis is that the resulting average profile is insensitive 
to the inherent asphericity and substructure (in projec- 
tion) of individual cluster mass distributions, as well as 
to uncorrelated large-scale structure projected along the 
same line of sight. Consequently, the statistical precision 
can be boosted by stacking together a number of clusters, 
especially on small angular scales (see Okabe et al. 2010). 

To do this, we first define a new set of radial bands 
in which the mass profiles of individual clusters are re- 
evaluated for a stacking analysis. Here we scale each 
cluster mass profile according to the cluster virial radius 
r v i r obtained by our non-parametric method (see § 4.5). 
For each cluster, we construct an M x N projection ma- 
trix Vji that projects the mass profile Ki (i = 1, 2, N) 
of the cluster onto the new radial bands scaled in units 
of r V i r (j = 1,2,..., M). Assuming a constant density in 
each radial band, the projection matrix Vji is uniquely 
specified by the conservation of mass. With this projec- 
tion matrix, the mass profile in the new basis is written 

15 Note, without the inner strong lensing information com- 
bined, Umetsu & Broadhurst (2008) found M vir = (1.38 ±0.14) X 
1O 15 M0 h _1 (the errors represent only the statistical one) from 
an entropy-regularized maximum-likelihood combination of Sub- 
aru distortion and magnification data sets, in excellent agreement 
with the one dimensional results in this work. 
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Figure 4. Non-parametric cumulative projected (blue-hatched) and spherical (gray-shaded) mass profiles, M2d(< 0) and M%y>(< 8), 
respectively, shown separately for the five clusters. For each cluster, the blue-hatched area in the upper panel shows the 68.3% confidence 
interval for M2d(< 6) at each radius estimated from a Monte-Carlo error analysis taking into account the error covariance matrix of the 
projected mass profile k(9) (see Figure 3). Similarly, the gray-shaded area shows the results for the deprojected mass profile M-jY>(< 8) 
assuming spherical symmetry. Shown in the lower panel is the mean spherical overdensity ~p~(< r) = A/3d(< r)/(47rr 3 /3) interior to radius 
r in units of the critical density of the universe p C rit(z), A(< r) = p(< r) / p CI i t (z), reconstructed from our comprehensive weak lensing 
analysis. The gray-shaded area shows the 68.3% confidence interval at each radius estimated from a Monte-Carlo error analysis taking into 
account the error covariance matrix of the projected mass profile k(8) (see Figure 3). The horizontal solid, dashed, and dotted lines show 
three representative values for the fractional overdensity, A = A v i r ,200, and 2500, respectively, for each cluster. See also Table 7. 



Table 7 

Three-dimensional Cluster Mass from a Non-parametric Deprojcction Analysis. 



Cluster 


A = 2500 

M2500 ^2500 

(IO^A/q/i- 1 ) (Mpch- 1 ) 


A = 500 

M500 r 5 oo 
(10 15 M r 1 ) (Mpch- 1 ) 


A = 200 

M200 V200 

(IO^Mq/i- 1 ) (Mpch- 1 ) 


A = A vir 

-^_/"vir T vir 

(lO^A/Qh- 1 ) (Mpcft- 1 ) 


A1689 
A1703 
A370 
C10024+17 
RXJ1347-11 


0.569 ±0.073 0.547 ±0.025 
0.404 ± 0.069 0.471 ± 0.028 
0.529 ±0.134 0.496 ±0.040 
0.462 ± 0.063 0.472 ± 0.023 
0.587 ±0.060 0.500 ±0.017 


0.883 ±0.117 1.083 ±0.050 
0.787 ±0.126 1.006 ±0.058 
1.315 ±0.161 1.152 ±0.046 
0.863 ±0.136 0.993 ±0.055 
0.949 ±0.145 1.003 ±0.051 


1.179 ±0.170 1.618 ±0.085 
1.143 ±0.196 1.546 ±0.097 
2.213 ±0.270 1.860 ±0.079 
1.194 ±0.188 1.502 ±0.081 
0.972 ±0.208 1.373 ±0.089 


1.300 ±0.205 2.011 ±0.113 
1.325 ±0.221 1.915 ±0.148 
2.399 ±0.249 2.215 ± 0.079 
1.329 ±0.224 1.799 ±0.105 
1.150 ±0.250 1.663 ±0.115 



Note. — M A 
luster rcdshift z j 



>(< »- A ) 



as 



k = Pk. 



Accordingly, the error covariance matrix in the new basis 
is 

C = VCV 1 . (21) 

With the mass profiles of individual clusters on a com- 
mon basis, we can stack the clusters to produce an aver- 
aged mass profile. Here we re-evaluate the mass profiles 
of the individual clusters in 8 logarithmically-spaced ra- 
dial bins over the range of radii R = [0.05, 1.7]r v i r . Since 
the noise in different clusters is uncorrclated, the mass 
profiles of individual clusters can be co-added according 
to (e.g., Sayers et al. 2009) 



where the index n runs over all of the clusters, and w n 
(20) ^ s t ne inverse critical surface mass density for the nth 
cluster, w n = (S-J t ) n . The error covariance matrix for 
the stacked mass profile (E) is obtained as 




(23) 



V n 



where the index n runs over all of the clusters. 

We show in the top panel of Figure 5 the re- 
sulting model-independent average mass profile (E) 
with its statistical la uncertainty as a function of 
the scaled projected radius i?/r vir , obtained by stack- 
ing the five clusters using equations (22) and (23) 
(see also Mandelbaum et al. 2006; Johnston et al. 2007; 
Sheldon et al. 2009). We note that the effect of different 
(22) cluster rcdshifts has been taken into account by proper 
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Figure 5. (Top): The model-independent average mass profile (T,)(R) (filled squares) with its statistical la uncertainty as a function of 
the projected radius R scaled with r v i r , which is obtained by stacking the lensing convergence profiles for the five clusters, A1689, A1703, 
A370, C10024+17, and RXJ1347-11, shown as solid lines. (Bottom): The logarithmic slope of the stacked mass profile (open squares with 
error bars), d ln{E)/d In R, is shown as a function of the scaled projected radius i?/r v j r , along with NFW model predictions with c v j r = 5 
and c v ; r = 10 for comparison. The projected logarithmic slope steepens from 720 = —1.01 ± 0.09 at R 0.1r v ; r to 72D = —1.92 ± 0.51 at 
R a 0.9r vir . 



error propagation in terms of the lensing efficiency func- 
tions (w n ) of individual clusters to average over. The 
stacked mass profile exhibits a fairly smooth radial trend, 
and is detected at a high significance level of 37a out to 
R = 1.7r vir . 

In the bottom panel of Figure 5, we plot the logarith- 
mic density slope 72d(-R) = dhx (£) /dlni? of the stacked 
mass profile as a function of the scaled projected radius 
along with NFW model predictions with c v i r = 5 and 
c v ir = 10. The logarithmic gradient of the average pro- 
file shows a slight steepening trend with increasing ra- 
dius in projection, consistent with NFW profiles with 
Cvir = 5 — 10. Finally, we quote model-independent con- 
straints on the average logarithmic density slope to be 
72D = -1.01±0.09atiJ« 0.1r vir and 7 2D = -1.92±0.51 
at R « 0.9r V i r for the five clusters. 

4.7. Combining Weak and Strong Lensing 

The Subaru data allow the weak lensing profiles of in- 
dividual clusters to be accurately measured in several in- 
dependent radial bins in the subcritical regime (9 > cin ). 
The projected mass profile can be unambiguously recov- 
ered on an individual cluster basis from combined weak- 
lcnsing shape distortion and magnification bias measure- 
ments. 

Here we combine our wcak-lcnsing profiles with de- 
tailed strong-lcnsing information for the inner ^ 200 kpc 
region of these clusters, for which we have identified 
many new sets of multiple images from deep HST/ ACS 
observations (Broadhurst et al. 2005b; Zitrin et al. 2009, 
2010), for a full determination of the entire mass pro- 
files of the five well-studied clusters. Figure 6 shows a 
sample of joint mass profiles for our five clusters recov- 
ered over two decades of radius ranging from 10 kpc h~ x 
to 3000 kpc ft, -1 . Note in this comparison we have ex- 
cluded the central weak-lensing bin k(< 9 m in) and the 
strong-lcnsing data points at radii overlapping with the 
Subaru data. In each case, the weak and strong lensing 
are in excellent agreement where the data overlap (typi- 
cally around R ~ 150 kpc /i -1 ), and the joint mass pro- 



files form well-defined radial profiles with a continuously- 
steepening radial trend from the central region to beyond 
the virial radius ( 1.7r v j r ). 

Our high-quality lensing data, covering the entire clus- 
ter, allow us to place useful constraints on the gNFW 
structure parameters (§ 4.3), namely, the central cusp 
slope a as well as the NFW virial mass and concentration 
parameters. Using our full lensing constraints, we obtain 
for each cluster the best-fit gNFW model as summarized 
in Table 6. For the halo mass and concentration param- 
eters (M v i r ,c_2), we find good agreement between the 
results with and without the inner strong-lcnsing profile 
combined (for C10024+17, see also Umetsu et al. 2010). 
Our joint mass profiles for the entire cluster region are 
consistent with a generalized form of the NFW density 
profile with modest variations in the central cusp slope 
(a ,$ 0.9), except for the ongoing merger RXJ1347-11 
(see Mason et al. 2010), for which we find an unaccept- 
able fit with a reduced \ 2 of 58 for 28 degrees of freedom, 
due to various local deviations from the model, especially 
in the inner mass profile which is tightly constrained by 
strong lensing. The best-fit value of a derived for the 
relaxed cluster A1703 is a = 0.93^0 257 being consis- 
tent with NFW (a = 1), which is in excellent agreement 
with independent lensing results by Oguri et al. (2009), 
a = 0.9±Q 4, and Richard et al. (2009), a = 0.92t^|. 

5. SUMMARY AND DISCUSSION 

Wc have developed a new method for a direct recon- 
struction of the projected mass profile of galaxy clusters 
from combined weak-lensing distortion and magnification 
measurements (§ 3.1 and § 3.2) within a Bayesian statis- 
tical framework, which allows for a full parameter-space 
extraction of the underlying signal. This method applies 
to the full range of radius outside the Einstein radius, 
where nonlinearity between the surface mass density and 
the observables extends to a radius of a few arcminutes, 
and recovers the absolute mass normalization. A proper 
Bayesian statistical analysis is essential to explore the 
entire parameter space and investigate the parameter de- 



14 



Bayesian Analysis of Weak Lensing Distortion and Magnification Measurements 



C\2 

o 

Oh 



Cluster mass profiles from ACS+Subaru lensing 



© 



10 



16 



,15 



,14 



W 1( 

•i-H 

ti 1, 

<u 1 ( 

X! 

m 
w 
cd 

a io 13 ^ 

cd 
o 
cd 



u 

CO 



10 



12 



10 



Strong lensing 



Weak lensing 




* A1689 (2=0.183) 

x A1703 (2=0.281) 

• A370 (2=0.375) 

★ C10024+1654 (2=0.395) 

■ RXJ1347-1145 (2=0.451) 



100 

R [kpc/h] 



1000 



- o 



'o 



- o 



Figure 6. Full mass profiles for A1689 (triangles), A1703 (crosses), A370 (circles), C10024+17 (stars), and RXJ1347-11 (squares) 
over a wide range of radius, from lOkpc/i -1 to 3000 kpc/i -1 , derived from ACS strong-lensing (R ?S 150 kpc /i _1 ) and Subaru weak- 
lensing (R ^ 150 kpc h" 1 ) measurements, showing a continuously steepening radial trend out to and beyond the cluster virial radius 
(r vir ~ 2Mpch~ 1 ). 



generacies (§ 3.3), arising from the mass-sheet degener- 
acy (§ 2). This method can be readily generalized for 
a statistical analysis using stacked lensing profiles of a 
sample of clusters. 

We have applied our comprehensive lensing method 
to a sample of five high-mass clusters for which 
detailed strong-lensing information is readily avail- 
able from EST I ACS observations. The deep Sub- 
aru multi-band photometry, in conjunction with our 
background-selection techniques (Mcdezinski et al. 2010; 
Umetsu et al. 2010), allows for a secure selection of un- 
contaminated blue and red background populations. In 
each cluster, a strong depletion of the red galaxy counts 
has been detected at a significance level of 6<7-11<t (Table 
5). A comparison shows clear consistency between two 
independent lensing observables with different system- 
atics, ensuring the reliability of our weak-lcnsing anal- 
ysis (Figure 1). The combination of independent Sub- 
aru distortion and magnification data breaks the mass- 
sheet degeneracy, as examined by our Bayesian statisti- 
cal analysis (Figure 2). Excluding magnification data, 
on the other hand, strongly modifies the Gaussian shape 
of the marginalized posterior PDFs, producing long non- 
Gaussian tails and broadening the distribution function, 
resulting in large reconstruction errors. The improve- 
ment here from adding the magnification measurements 
is significant, ~ 30% in terms of cluster mass profile mea- 
surements (Table 5). 

We have formed a model-independent mass profile 
from stacking the clusters, which is detected at 37tr out to 
beyond the virial radius, R sa 1.7r v ; r . We found that the 
projected logarithmic slope, 72d(-R) = din E(i?)/dlni?, 
steepens from 72D = —1.01 ± 0.09 at R « 0.1r v j r to 
72D = —1-92 ± 0.51 at R « 0.9r v i r , consistent with 



NFW profiles with c v - n - = 5 — 10. We also obtained 
for each cluster inner strong-lensing based mass profiles 
from deep HST/ ACS observations, which we have shown 
overlap well with the outer Subaru-based profiles and to- 
gether are well described by a generalized Navarro- Frenk- 
White profile, except for the ongoing merger RXJ1347-11 
(Mason et al. 2010), with modest variations in the cen- 
tral cusp slope (a ^ 0.9), perhaps related to the dynam- 
ical state of the cluster. 

These high-mass lensing clusters with large Einstein 
radii appear to be centrally concentrated in projection, as 
found in several other well studied massive clusters from 
careful lensing work (Gavazzi et al. 2003; Kneib et al. 
2003; Broadhurst et al. 2005a; Limousin et al. 2007; 
Broadhurst et al. 2008; Broadhurst & Barkana 2008; 
Oguri et al. 2009; Umetsu et al. 2010; Zitrin et al. 2011). 
An accurate characterization of the observed sample is 
crucial for any cluster-based cosmological tests. It has 
been suggested that clusters selected with giant arcs 
represent a highly biased population in the context of 
ACDM. For those clusters with large Einstein radii (say, 
t9 c i n > 20"), a large statistical bias of about (40 — 60)% is 
derived from iV-body simulations (Hennawi et al. 2007; 
Oguri & Blandford 2009), representing the most triax- 
ial cases. The mean level of mass concentration inferred 
from our weak lensing analysis is high, (c V j r ) = 7.9 ± 0.8 
(simply ignoring the mass and redshift dependencies), for 
our sample with (M vir ) = (1.48 ± 0.10) x 1O 15 M h' 1 , 
as compared to the ACDM predictions, (c v j r ) = 3.4^'q 
(the errors quoted represent a la lognormal scatter of 
fj(log 10 Cvir) = 0.15; see Duffy et al. 2008), evaluated at 
the median redshift of our sample, Zd — 0.38. This rep- 
resents an overall discrepancy of 5cr with respect to the 
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predictions, without taking into account the effects of 
projection bias. This apparent discrepancy is also evi- 
dent when the weak-lensing mass profiles are combined 
with the inner strong-lcnsing information from deep ACS 
observations (§ 4.7). Applying a bias correction of 50% 
(see Broadhurst et al. 2008; Oguri et al. 2009), the dis- 
crepancy in (c v ir) is reduced to the 3tr level. 

Another possible source of systematic errors in the 
mass profile determination is the cluster off-centering 
effect (§ 4.3). The effect is essentially the smoothing 
of the central lensing signal (see Oguri & Takada 2011), 
which flattens the recovered convergence profile below 
the offset scale, and therefore reduces the derived mass 
concentration (c_ 2 ) and cusp slope (a) parameters. For 
our sample of clusters, we found a cluster centering off- 
set of typically A9 ^ 5" using our detailed strong-lcnsing 
models. This level of centering offset (Ai? ^ 20kpc/i _1 ) 
is much smaller than the estimated values for the inner 
characteristic radius of our clusters, r_2 ~ 200kpc/i~ 1 , 
and hence may not significantly affect our concentration 
measurements. However, this could potentially lead to 
an underestimation of the central cusp slope a. 

Our results are consistent with previous lensing work 
which similarly detected a concentration excess in the 
lensing-based measurements for strong-lensing clusters 
(Comerford & Natarajan 2007; Broadhurst et al. 2008; 



Oguri et al. 2009). Our findings could imply, for these 
clusters, either substantial mass projected along the line 
of sight, due in part to intrinsic halo triaxiality, or a 
highcr-than-expected concentration of dark matter. Nev- 
ertheless, the overall level of uncertainties may be too 
large to robustly test the CDM predictions with the 
present sample. The forthcoming space telescope clus- 
ter survey, CLASH 16 , will provide a definitive deriva- 
tion of mass profiles for a larger, lensing-unbiased sam- 
ple of relaxed X-ray clusters (> 5kcV), combining high- 
resolution panchromatic space imaging with deep, wide- 
field Subaru weak-lensing observations, for a definitive 
determination of the representative mass profile of mas- 
sive clusters. 
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APPENDIX 
A. MASKING CORRECTION 

For the cluster magnification analysis, the masking effect due to bright cluster galaxies, bright foreground ob- 
jects, and saturated objects has to be properly taken into account and corrected for (Broadhurst et al. 2005a; 
Umetsu & Broadhurst 2008; Umetsu et al. 2010). Here we describe two different approaches for estimating the ef- 
fect of masking in the magnification bias measurement. 

A.l. Method A 

In this method (Method A), we describe each masking object as an ellipse specified by its structure parameters, such 
as the size, axis ratio, and position angle. SExtractor provides such useful estimates of object properties (magnitude, 
size, axis ratio, position angle, and so on) for calculating masking areas and the resulting correction. In practice, 
we account for the masking of observed sky by excluding a generous area nab around each masking object, where a 
and b arc defined as ^ mas k(= 2 — 4) times the major (A_IMAGE) and minor axes (B_IMAGE) computed from SExtractor, 
corresponding roughly to the isophotal detection limit (see Umetsu & Broadhurst (2008) and Umetsu et al. (2010)). 
We then calculate the correction factor for this masking effect as a function of radius from the cluster center, and 
renormalize the number density accordingly. Note that the actual choice of ^ mas k is insensitive to the resulting 
magnification measurement n^/no = /i 2 5s_1 (Umetsu & Broadhurst 2008; Umetsu et al. 2010), as this effect will 
cancel out when taking a ratio of the local surface density to the unlensed mean number density no (see § 5.5.3 of 
Umetsu & Broadhurst 2008). However, this method is not suitable for estimating masking areas by saturated stars and 
stellar trails, which cannot be described by simple ellipses. For this case, we manually place a rectangular region mask 
around each saturated object to exclude a proper area. This approach has been commonly adopted in previous studies 
of cluster magnification bias (Taylor et al. 1998; Broadhurst et al. 2005a; Umetsu & Broadhurst 2008; Umetsu et al. 
2010). 

In the top panel of Figure 7, we show as an example the masked region in A1689 based on the Subaru i'-band 
photometry, as obtained using Method A (see also Umetsu <fc Broadhurst 2008). The gray area shows bright masking 
objects (i' < 22 AB mag), each described as an ellipse. The black rectangular regions are masked areas for saturated 
stars and bright stellar trails. We show in Figure 8 the level of mask area correction in A1689 as a function of radius 
from the cluster center. The masking area is negligible (a few %) at large radii, and increases up to ~ 20% of the sky 
close to the cluster center. Typically, the level of masking correction is found to be sufficiently small (^ 1), so that 
the uncertainty in the masking correction is of the second order, and hence negligible. 

A.2. Method B 

Here we describe an alternative method (Method B) for the masking correction which we have developed in this 
work. This method has practical advantages in comparison to Method A: that is, it is easy to implement and can be 



16 Cluster Lensing And Supernova Survey with Hubble (P.I.: M. 
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Figure 7. Distribution of masked regions in A1689 (z = 0.183). Superposed are the concentric radial bands used for calculating the radial 
profiles centered on the BCG. Top: Results obtained using Method A (sec Appendix A.l). The gray area shows bright masking objects 
(£' < 22 AB mag), each described as an ellipse. The black rectangular regions are masked areas for saturated stars and bright stellar trails. 
Bottom: Results obtained using Method B (see Appendix A.2) based on SExtractor's check-image output with CHECKIMAGE.TYPE = OBJECT. 
The black area shows those connected pixels that belong to bright foreground objects, bright cluster galaxies, saturated stars, and stellar 
trails. 
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Figure 8. The azimuthally-avcragcd correction factor / mas k(0) for the masking effect as a function of radius from the center of A1689. 
The red dashed curve shows the level of correction required due to masking, as estimated by Umetsu & Broadhurst (2008) using Method 
A (Appendix A.l). The black solid curve shows the results using Method B (Appendix A.2). The two methods give consistent results in 
terms of both shape and amplitude of the radial profile of the correction factor. 

fully automated once the configuration of analysis parameters is tuned, especially useful for large scale sky surveys 
such as the Subaru HSC survey. 
In this approach, we use a pixelized flag image to exclude masking regions, namely those connected pixels that 
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belong to bright foreground objects, bright cluster galaxies, saturated stars, and stellar trails. This can be done using 
SExtractor's check-image output with CHECKIMAGE_TYPE = OBJECTS. We tunc SExtractor configuration's detection 
parameters in order to optimize the detection for a particular object size and brightness. In our analysis of Subaru 
images, the per-pixel detection threshold above the local sky background is set to 5a (DETECT_THRESH = 5), and 
the minimum number of connected pixels above threshold is set to 300 (DETECT_MINAREA = 300) with 0.202" pixel -1 
sampling. The bottom panel of Figure 7 shows the distribution of masked regions in A1689 as obtained using Method 
B. Figure 8 compares the masking correction factor in A1689 obtained using Methods A and B, shown as a function 
of radius from the cluster center. It is shown that the two methods give consistent results in terms of both shape and 
amplitude of the radial profile of the correction factor. 

B. DISCRETIZED ESTIMATOR FOR THE AVERAGED CONVERGENCE 

In this Appendix, we aim to derive a discrete expression for the mean interior convergence 7t(< 9) using the 
azimuthally-averaged convergence n(6). In the continuous limit, the mean convergence k(< 9) interior to radius 8 
can be expressed in terms of k(9) as 



R(0) = J- / d\n9'8' 2 K{9'). 
v Jo 



(Bl) 



For a given set of (N + 1) annular radii 9[ (I = 1,2, N + 1), defining N radial bands in the range 6 m i n = Q\ < 6 < 
9n+i = #max, a discretized estimator for k(< 9) can be written in the following way: 



«(< °i) = %^(< 0«nm) + Aln^«(«,), (B2) 

with Aln#; = (9i + i — 0i)/9i and 9i being the area- weighted center of the zth annulus defined by 9i and in the 

continuous limit, we have 



_ 2 9f + 6 i+1 + 9i9 i+ i 
3 9i + 9i+i 

C. DISCRETIZED ESTIMATORS FOR THE LENSING PROFILES 



(B3) 



We derive expressions for the binned tangential distortion g + and magnification \i in terms of the binned convergence 
k, using the following relations: 

g+ m= ~ << /' , -f\ (CD 

1 - K{9i) 

p(Si) = t , (C2) 

(i-«(fl*)) a (i-fl3.(fl0) 

where both the quantities depend on the mean convergence k interior to the radius 9i, which is the center of the ith. 
radial band of [9i, 9i + {\ (see Appendix B). By assuming a constant density in each radial band and by noting that 9i 
is the median radius of the zth radial band, k(< 9i) can be well approximated by 

«(<3O = 5(«(<0O + «(<0t+i)) (C3) 

where k(< 9i) and k(< 9i+i) can be computed using the formulae given in Appendices B and C. 

If we are to combine background samples of different populations (a and b) and different lensing depths (/3), we 
should take into account population-to-population variations in the lensing depth, u>i, a = (/3}b/(/3) a , which can be 
measured from deep multi-band imaging, such as the 30-band COSMOS database (Ilbert et al. 2009). The tangential 

distortion g^t' for Population (b) can be expressed in terms of the convergence for Population (a) as 

3+ - ■ (C4) 

Similarly, the magnification jjf^ for Population (b) can be expressed in terms of u}b a tt^ . 
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d. supplemental material 

In this supplemental material section, we compare observed tangential-distortion (shear) g+ and count-depletion 
(magnification) radial profiles with their respective reconstructed profiles, taking A370 at z c i = 0.375 (Figure 9) 
and C10024+17 at Zd = 0.395 (Figure 10) as examples. We note that this supplemental material is added for the arXiv 
version only, and is not included in the published ApJ version (Umetsu et al. 2011, ApJ, 729, 127). 
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Figure 9. Left: wcak-lcnsing radial profiles of A370 [zd = 0.375) as measured from background galaxies registered in Subaru BjR c z' 
images (Table 1). The top panel shows the tangential reduced shear profile g+(0) (squares) based on Subaru distortion data of the full 
background (red+blue) sample (Table 3). The bottom panel shows the count depletion profiles n^(9) due to magnification for a flux-limited 
sample of red background galaxies. The circles and triangles show the respective results with and without the mask correction due to 
bright foreground objects and cluster members. The horizontal bar represents the constraints on the unlensed count normalization, no, as 
estimated from Subaru data. Also shown in each panel is the joint Bayesian fit (gray area, 68% CL) to both profiles, Right: the joint mass 
profile solution of the shear and magnification data sets shown in the left panels, obtained following the Bayesian method described in the 
main text (Section 3), effectively breaking the mass-sheet degeneracy. The top panel shows the reconstructed surface mass density profile 
k = S/S cr i t (squares), in good agreement with the standard NFW form (gray area, 68% CL; see Table 6). The innermost bin represents 
the average surface mass density re m i n = re(< # m i n ) interior to the inner radial boundary f3 m i n = 0.7' of the wcak-lcnsing distortion data. 
The bottom panel shows the projected cumulative mass profile M20 (red), along with the best-fit NFW model (gray area, 68% CL; Table 
6). 
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Figure 10. Same as in Figure 9, but shown for the cluster C10024+17 (z,i = 0.395). A strong depletion of the red galaxy counts is clearly 
detected (bottom left panel) down to the inner radial boundary of the weak lensing data, f9 m i n =0.6'. A consistent mass profile solution 
s = {75 m ; n , Ki} i=1 (top right panel) is obtained from a joint Bayesian fit to both lensing profiles shown in the left panels. 
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